Thermally- Activated Magnetic Reversal Induced by a Spin-Polarized Current 
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We have measured the statistical properties of magnetic reversal in nanomagnets driven by a spin- 
polarized current. Like reversal induced by a magnetic field, spin-transfer-driven reversal near room 
temperature exhibits the properties of thermally-activated escape over an effective barrier. However, 
the spin-transfer effect produces qualitatively different behaviors than an applied magnetic field. We 
discuss an effective current vs. field phase diagram. If the current and field are tuned so that their 
effects oppose one another, the magnet can be driven superparamagnetic. 

PACS numbers: 73.40.-c, 75.60.Jk, 75.70.Pa 
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A spin-polarized current traversing a magnetic multi- 
layer can, through exchange interactions, alter the ori- 
entation of ferromagnetic moments, producing domain 
reversal or exciting spin waves. Following early predic- 
tions [0-^) , these spin-transfer effects have been observed 
j| and they are generating interest as an alternative 
to the use of magnetic fields for switching elements in 
magnetic memories. Competing theories of the effect 
have thus far considered only the limit when the tem- 
perature T = 0. The initial theories, which we shall 
call the torque model JlPJl^], predict that domain re- 
versal should occur when spin transfer from the current 
produces a torque that exceeds the magnetic damping. 
The switching current is therefore not determined by an 
energy barrier. Alternative approaches have predicted 
that the spin-polarized current may provide an effective 
magnetic field favoring parallel or antiparallel alignment 
of adjacent magnetic layers, thereby altering an energy 
barrier for reversal [p~3|JT^ ] . Here we report, for mea- 
surements of spin-transfer-driven reversal near room tem- 
perature, broad distributions of switching currents that 
depend strongly on temperature, similar to the famil- 
iar distributions of switching fields measured when an 
applied magnetic field drives thermally-activated mag- 
netic reversal ]l5|JT^ |. This indicates that the primary 
effect of spin transfer near room temperature is to al- 
ter a thermally-activated over-barrier switching process. 
Nevertheless, by comparing the effects of applied cur- 
rents and applied magnetic fields, we demonstrate that 
the spin-transfer effect cannot be understood as an effec- 
tive magnetic field favoring parallel or antiparallel mag- 
netic alignment pl[l7j]. If we extrapolate our switching 
currents and fields to T — 0, we find qualitative agreement 
with the torque model. 

Figure 1 (a) is a schematic of our device geometry. The 
fabrication process Q employs electron beam lithogra- 
phy and ion milling to form a pillar with cross section 
ranging from ~ 50 x 50 nm (sample 1) to ~ 130 x 60 nm 
(sample 3) from a multilayer of 80 nm Cu/40 nm Co/6 
nm Cu/3 nm Co/10 nm Au. The milling step is timed so 
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FIG. 1. (a) Schematic of the nanopillar device, (b in- 
set) Differential resistance vs. / for device I. H is along the 
easy axis of the nanomagnet. (b) Distribution of switching 
currents for device 1, for the parallel to antiparallel transi- 
tion, (c) Waiting time distributions for device 2. The dis- 
tributions are fit to the function e~ l ^ T . (d) Dependence of 
the mean switching current (I c ) on T (closed squares) and 
current-sweep rate at room T (open squares) for device 3 at 
H = 0. The bars show the width 2<7r c of the distributions. 

that the thicker Co layer is left as an extended film. The 
differential resistance dV/ dl as a function of bias current 
/ perpendicular to the layers is plotted in the inset to 
Fig. 1(b). If the moments of the two magnetic layers are 
initially parallel (P) and the current is swept from neg- 
ative to positive | |Djl , dV/dl jumps at a critical current 
Ij, where the moment of the small Co nanomagnet is 
driven antiparallel (AP) to the thicker Co film. The de- 
vice remains in this AP state until the current is swept 
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down past a negative value I~ , at which point the spin- 
transfer effect drives the nanomagnet back parallel to the 
thicker film. This asymmetry in / is in agreement with 
the theories HMFlM| ,|9[-|j] , and with previous switch- 
ing studies |||8|j£ ^|Tl ] . 

The value of current at which the magnet reverses 
varies from sweep to sweep. A histogram of at room T 
is shown in Fig. 1(b) for a current-sweep rate of 80 fiA/s. 
Similar histograms are found for all of the eight samples 
we have studied in detail. The stochastic nature of the 
switching is confirmed in experiments in which we hold 
the sample at a fixed I near a switching threshold; there 
is a waiting time before switching that also displays a 
broad distribution |H. Probabilities P(t) that the mag- 
net has not reversed in a time t are plotted in Fig. 1(c) for 
a second sample and compared to exponential decay. As 
a function of either decreasing T or increasing sweep rate 
dl/dt, the distributions of I c shift strongly to larger val- 
ues of |/| (Fig. 1(d)). At temperatures below 100-150 K, 
the switching events generally consist of multiple jumps 
in resistance rather than a single jump. Related multiple 
jumps have been observed previously [||||, and they indi- 
cate that the reversal mechanism is not coherent rotation 
but rather a more complicated process like domain wall 
nucleation 
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Our main focus in this Letter will be on the statistical 
distributions of switching currents and switching fields 
as both / and magnetic field H are applied to the sam- 
ple. By studying the interplay between spin transfer and 
if, we can gain new insights into the mechanism under- 
lying spin transfer. Fig. 2(a) shows the mean values of 
the current-induced reversal distribution as a function of 
constant H at room T measured with a current-sweep 
rate of 80/zA/s. H for all data in the paper is applied 
within 5° of the easy axis of the nanomagnet. Both criti- 
cal currents, and shift towards more positive / 
as a function of H, with the AP— »P transition even- 
tually shifting more rapidly until it intersects (I^f) an d 
the nanomagnet is no longer bistable ||. In Fig. 2(b) 
we show the corresponding standard deviations for the 
switching currents, af and aj . These show unexpected 
differences. At low fields, <t/~ and a J are roughly equal 
to one another. However, as H increases near the values 
required for a field-driven transition to the parallel state, 
<jj increases by more than threefold, while af decreases 
slightly. 

If the roles of I and H are reversed, the results give 
a first indication that spin transfer acts quite differ- 
ently than the magnetic field. The inset to Fig. 2(c) 
shows a plot of dV/dl near zero current bias versus H 
along the easy axis for the same device as that used 
for Figs. 2(a) and 2(b). The low-field transition to the 
higher-resistance state is due to the reversal of the lower- 
coercivity extended Co film, while the higher-field tran- 
sition to low resistance corresponds to the switching field 



(a) 

< 
E 

A 
V 



0.5 



0.0 



-0.5 



(b) 



0.02 



< 
£ 



0.01 



0.00 



(C) 150 
P 

S 100 



□ 


□ □ D 


□ 


Pto AP 


• 


AP to P 


• 


• • • 



□ □□□□□□□□□ nyC 



T = 293 K 



■* m B - 




□ 


Pto AP 


• 


AP to P 



20 40 60 80 100 120 140 

(mT) 



v 



50 



+ 
i \ 



T = 293 K 



H.H (mT) ~ ■ 



> 

N- 

"+-+. 



"+ -+ 



2.5 
2.0 

P 
1.5 £ 

1.0* 



-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

l(mA) 

FIG. 2. (a) Mean and (b) standard deviation of I c for 
device 1 as a function of H . Lines in (b) are fits of the stan- 
dard deviations to Eq. (4). (c) Mean and standard deviation 
of H c as a function of / for device 1. Inset: the low-bias 
magnetoresistance loop. 

of the nanomagnet, H c . (H c ) and cth c as a function of 
/ are shown in Fig. 2(c). While (H c ) monotonically in- 
creases with /, there is a distinct peak in an c , occurring 
at -0.2 mA. 

The striking behavior exhibited in Fig. 2 can be ex- 
plained by extending the theory of thermally-activated 
magnetic-field-induced switching. We will employ the 
Kurkijarvi approach, in which the effective activation 
barrier to reversal is tuned by varying an external param- 
eter at a constant rate [Q. For generality we will call 
the barrier-reducing parameter U, and we will apply the 
theory when D = I as well as D = H . The Kurkijarvi cal- 
culation assumes that the activation barrier has the ap- 
proximate form of a power law U (D) = C/ (1 — D /D c o) aD 
where D c q is the parameter at which the barrier vanishes 
at T = 0, and et£> is a constant. The mean and standard 
deviation of the distributions for the switching point D c 
are, to leading order pa, 
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where D is the sweep rate and To is the attempt time. 

We should note that, in addition to reducing the bar- 
rier to reversal, the application of current could also heat 
the device. However, by comparing the magnitude of I- 
dependent changes in DC resistance at the relevant cur- 
rent levels (a few mA) to the T-dependence of the low- 
bias resistance, we estimate that the devices are heated 
by at most 2-3 K above room temperature. This corre- 
sponds to a 1% effect on (J c ), which can be neglected in 
our analysis. 

If we apply a fixed magnetic field while sweeping I, 
H has two effects on the I c distributions. First, H al- 
ters the effective zero-current barrier height U$. For 
the AP— »P transition, the barrier is reduced from Uq 
to U Q {1 - H/H c0 ) aH , while the P->AP barrier is in- 
creased to Uq(1 + H / H c o) aH . Second, the magnetic field 
can modify the zero-temperature critical currents 1^ for 
the two transitions; the form will depend on the mi- 
croscopic model. However, in the models proposed to 
date Jl|,pp^], this dependence is linear, so that we will 
take I%{H) = I^(0)(1 - H/Hf), where Hf are model- 
dependent. Inserting these two quantities into Eq. (2) 
and neglecting the weak ii-dependence of A yields 



a%{H) k (1 - H/Hf)/{1 ± H/H c0 ) 



a H /ai 



(4) 



If spin transfer merely acted as an additional effective 
field in the direction of H, then we should have aj = an 
and Hf 1 = TH c o- In this case the numerator and denomi- 
nator in Eq. (4) cancel and af c should be f/-independent. 
This does not describe the data. In contrast, within 
Slonczewski's torque model, H c q and \H^\ differ. For 
a thin- film nanomagnet, H c o is set by the small in-plane 
anisotropy H coerc i ve w 150 mT, while the field intercept 
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Hdemag, where H demag represents the 
additional effect of the demagnetizing field as the mo- 
ment precesses out-of-plane (for single-domain magnets 
undergoing coherent rotation, ^oHdemag — « 850 

mT for Co) |s|j2"6[]. As a result, within this approach a J 
diverges at H — H c0 , while aj~ slowly decreases, in ex- 
cellent agreement with the data. The lines in Fig. 2(c) 
illustrate the results of this model using /j,qH c o = 150 
mT, \[ioHg\ = 230 mT, aj = a#, and the scale factors 
af o (H = 0) = 0.007 mA, aJ {H = 0) = 0.0065 mA. The 
fact that \Hg\ is less than M/2 may be due to non-single- 
domain dynamics. 

The dependence of cth c on / can be understood within 
the same model by considering the nature of the T = 
stability boundaries for P and AP alignment. A sim- 
ple relation allows us to connect the measured histogram 
means and widths to the T — stability boundaries: from 
Eqs. (l)-(3), since the function A depends weakly on its 
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FIG. 3. (a) T = current-field stability diagram (crosses) 
for device 1, showing the parallel (P), antiparallel (AP), and 
bistable (P/AP) regimes, determined as described in the text. 
The circles and squares are the same data as in Fig. 2, and 
solid lines are guides to the eye. (b) Solid lines: Form of the 
stability diagram in the Slonczewski torque model for a sin- 
gle domain magnet with easy-plane and in-plane anisotropies. 
NSS is a regime where non-static states exist. The critical cur- 
rents 7^ are related through I^ /I~ = g(n)/g(0), where g is 
the Slonczewski polarization factor [2] . Dotted lines: Stability 
boundaries in the effective-field model. 

variables, to a good approximation af oc |J c o — (1^)1 an d 
uh c oc H c q — (H c ). We can estimate the proportional- 
ity constant self-consistently as follows. The normalized 



sweep rate 



0.05-0.1 s 1 and we use an attempt 



time To w 100 ns |W|,|27|]. By fitting the dependence of 
(7 C ) on sweep rate for this sample to Eq. (l)-(3), assum- 
ing the barrier exponents to be a/ = an ~ 1.5 pq] , we 
find an H = effective barrier of Uq — 1.5 - 2 eV. Insert- 
ing these values into Eq. (l)-(3) yields a proportionality 
constant ~ 0.1, with the dominant uncertainty associated 
with to- The resulting estimates for the T = stability 
boundaries, extrapolated from the room-T measurements 
of (/*) and (H c ), are marked by the crosses in Fig. 3(a). 
The maximum in ajj c is associated with the knee in the 
stability diagram where the critical-current line I~ (H) 
joins the critical- field line H c q(I). In this region the mag- 
net is maximally subject to thermally-activated reversal, 
through the combined effects of I and H, and therefore 
<jh (I) is a maximum. If spin transfer worked as an ef- 
fective field, Icq(H) and H c q{I) would fall on one line, 
and there would be no maximum in uh c {I) ■ 

In order to compare these results to the torque model, 
we have calculated the stability boundaries within the 
T = Slonczewski picture by numerically integrating 
the Landau-Lifschitz-Gilbert equation with a spin-torque 
term for a single-domain magnet with easy-plane and in- 
plane anisotropies (Fig. 3(b)). Although we do not ex- 
pect this model to be quantitatively accurate if the rever- 
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shape of the T = stability diagram for P and AP ori- 
entations, and (3) the distinct difference between the ef- 
fects of / and H on switching times in the telegraph- noise 
regime. Our data do not rule out a small effective-field 
contribution in addition to the torque term (l2|,[l4j ■ 
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FIG. 4. Inset: DC resistance versus time for device 1 with 
HoH = 140 mT and I = +0.51 mA at room T. (a,b) Mean 
switching times for the two resistance states. The time scales 
in (a) and (b) do not match precisely because of a small shift 
in the sample between measurements. 

sal mechanism is not single-domain, it has the qualitative 
features needed to understand the data. It predicts that 
the I~ (H) and H c q(I) lines are distinct, and intersect at 
a knee located at a negative value of I. 

If we bias the sample near the point /iqH = 140 mT 
and I — +0.5 mA where the spin-transfer effect and 
H oppose each other, we observe telegraph-noise-type 
switching between resistance states (inset, Fig. 4(a)) [Q. 
Unlike previous telegraph-noise studies in nanomagnets 
[ p5| , which were done by applying H perpendicular to 
the easy axis so that the moment jumped between two 
closely-separated angles, the jump here is between ap- 
proximately full P and AP alignment. Most remarkably, 
the mean switching times for the two types of transi- 
tions depend very differently on H and /. If H is held 
fixed along the easy axis and I is increased (Fig. 4(a)), 
tp^ap decreases exponentially, while tap^p increases 
only slightly. Varying H while holding / fixed (Fig. 4(b)), 
on the other hand, decreases tap^p exponentially, while 
t p—>ap increases much more slowly. These differences 
provide independent evidence that the spin-transfer ef- 
fect and H alter the switching of the nanomagnet in 
largely independent ways. By tuning / and H along the 
easy axis, both switching times can be shortened until 
the nanomagnet is effectively superparamagnetic. 

In summary, magnetic reversal driven by spin- 
polarized currents exhibits statistical properties of ther- 
mal activation over an energy barrier. This might appear 
to favor an effective-magnetic-field model [ p^JT^ ] over the 
torque model [QJ^]. However, our data show that the 
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torque model provide natural explanations for (1) the 
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